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Abstract 

5_i ■ Processes of the form pp — > anything — > X{Xj — > xx + yy (+$) are studied 

via a technique that may be viewed as an adaptation of time-honoured Dalitz 
plot analyses. Xi and Xj are new heavy states (with i,j = l...n), which may 
be identical or distinct; and xx and yy are necessarily distinct Standard Model 
(SM) fermion pairs whose invariant masses can be measured. A Dalitz-like 
plot of said invariant masses, M(xx) vs. M(yy), exhibits a topology connected 
to the masses and specific decay chains of Xi and Xj. Aside from relatively 
minor details, observed patterns consist of a collection of box and wedge shapes. 
This collection is model-dependent: comparison of the observed pattern to the 
possibilities for a specific model yields information on which new particle pair 
combinations are actually being produced, information beyond that extractable 
from conventional one-dimensional invariant mass distributions. The technique 
is illustrated via application to the Minimal Supersymmetric Standard Model 
(MSSM) process pp —>■ gg,gq,qq — > xfiXj e + e~ + fJ> + u~ {+$)■ Here the 
heavy states are neutralinos Xi (* = 2, 3, 4) — note Xi is excluded — which 
are produced in gluino/squark (g/q) cascade decay chains. Even with fairly 
modest expectations for the LHC performance during the first few years, this 
method still provides substantial insight into the neutralino mass spectrum and 
couplings if gluino/squark masses are relatively low (~ 400 GeV). 



*present address: Department of Physics, Sichuan University 



Introduction 



The field of particle physics is nearing a critical juncture: up to now the highly 
successful SM — whose predictions of various cross-sections and precision observables 
are in excellent agreement with data from the most advanced particle accelerators 
to date — has been sufficient to meet the experimental demands; however, the SM 
is theoretically incomplete and cannot continue to describe physics at energies much 
higher than 1 TeV. Most theoretical extensions of the SM designed to address this 
problem predict new heavy degrees of freedom at or near the TeV-scale. The soon-to- 
be-completed LHC, with a centre-of-mass energy of 14 TeV, should readily produce 
such heavy particles if they couple significantly to SM ones. Then experiment will 
certainly require more guidance than the SM can provide. Different extensions to the 
SM differ in the predicted number and types of new heavy particles. It is therefore 
imperative to understand what the decays of such heavy particles (at least some 
of which are typically unstable) would look like at a hadron collider in as model- 
independent a way as possible. 

This work is particularly concerned with neutral heavy particles produced in 
pairs, XiXj, with i, j = l...n — the exact value of n being model dependent, but for 
the present work it may be any integer greater than 1 (thus Xi and Xj may or may 
not be distinct). These pairs may be produced directly as per pp — > X^Xj or as a 
result of cascade decays from the production of other even heavier new particles (in 
fact the latter production mode is dominant in the specific case examined below). 

The introduction of new heavy particle states often comes with the introduction 
of a new conserved quantum number (or numbers) associated with a new discrete Z n 
symmetry (or symmetries) — for example, the attractive .R-parity |Tj conservation 
in numerous SUSY extensions 1 to the SM. Another example is found in some little 
Higgs models in which an extra symmetry is introduced to tame excessive flavour- 
changing neutral processes [2|. Conservation of such a new quantum number (s) 
typically dictates the pair production of new particle states as well as the stability of 
the lightest new particle which is "odd" under the new symmetry (features typically 
associated with i?-parity conserving SUSY scenarios but in fact they seem to be more 
generally applicable). 

Any sample of events collected over time may be a superposition of different X^Xj 
channels. The technique introduced here is ideally suited for precisely this situation. 
Unlike at an e + e~ collider, at a hadron collider the centre-of-mass energy of the 
parton-level hard scattering process cannot be controlled, and thus said parton-level 
centre-of-mass energy cannot be incrementally raised to scan through the different 
XiXj thresholds. Rather, all such channels may be produced simultaneously and 
must subsequently be disentangled to the extent it is possible in the decay analysis. 

The new heavy particles are assumed to decay (possibly indirectly) into pairs of 
SM fermions Xi — > yy (accompanied in some models containing stable but unde- 
tectable heavy particle states by the observation of substantial missing energy in the 
detector). Thus pair production and subsequent decay of the new heavy particles 

henceforth the acronym 'SUSY' will be used for both 'supersymmetry' and 'supersymmetric'. 



can result in end-states of the form xx + yy , where xx and yy are distinct SM 
particle pairs whose invariant masses are measurable with sufficient precision. For 
example, same-flavour oppositely-charged lepton pairs e + e" and (utilised in 

the following application to the MSSM [I]) might be chosen since these are most eas- 
ily extracted from the overwhelming QCD backgrounds at a hadron collider. Other 
choices, such as bb or t~t + , are also possible though, and might prove more appro- 
priate in some cases. Decays to SM gauge bosons may merit attention, though with 
decays to Z 0, s a Z°-veto to reduce backgrounds is no longer possible while decays 
to W's will require reconstruction of hadronically-decaying W's. The remainder of 
this work concentrates on decays into pairs of fermions, and, more specifically, into 
electrons and muons. Use of similar ideas for the pair-production of charged states 
(X^Xj) also might merit future investigation. 

Speaking generally, the experimentally measurable quantities of interest are the 
fermion pair invariant masses M(xx) and M(yy). Other processes besides the sought- 
after heavy particle decays may also produce an (xx, yf/)-topology. Thus cuts will 
probably be needed to purify the event sample, and a partially-contaminated sample 
may have to suffice. It will be shown below that, for several realistic MSSM scenarios 
including both signals and backgrounds, making a two-dimensional Dalitz-like 3 plot 
of M(xx) vs. M(yy) can reveal information about the spectrum of the heavy particles 
produced (kinematics) as well as relative production cross sections (dynamics). 

Topological Analysis 

Any Dalitz-like plot of M(xx) vs. M(yy) resulting from heavy particle pair produc- 
tion will be a superposition of specific topological shapes. At the coarsest level, these 
shapes may be bifurcated into two types: 

• Box-like — A 'box' in the M(xx)-M(yy) plane results from the decay 



since the invariant masses M(xx) and M(yy) are bounded from below by the 
masses of x, y (approximately zero if x, y are leptons) and above by the maxi- 
mum for M Xi (this is a well-defined limit if the model in question completely 
accounts for f, by particles which do not decay in the detector). 

• Wedge-like — A 'wedge' or 'L-shape' results from the decay 



i.e., if the xx came from X{ then the yy presumably comes from Xj — here 
it is assumed neither x- nor y-flavour number is violated in the heavy-particle 
decays 2 . Therefore < M(xx) < M x% - $ and < M(yy) < M x . - $. On 
the other hand if the decays are swapped then < M(xx) < M Xj — $ and 
< M(yy) < M Xi — $; the superposition of these two strips forms the wedge. 

2 It is possible, for example, to have e ± fi zfl decays from neutralinos in the lepton-flavour-conserving 
MSSM, but the branching ratios (BRs) for such decay modes are small and generally negligible. 



XiXi ^xx + yy {+$) , 



(1) 



XiX s ->xx + yy {+$) (i ^ j) ; 



(2) 



The manner in which Xj and Xj decay may introduce new features on top of 
these two basic forms. For example, whether the decay proceeds through a series of 
two-body decays or via a three-body decay Furthermore, if some Xi involved in the 
decay chain has two or more ways to decay to xx and yy; e.g., if two or more decay 
chains resulting in Xi — ■> Xj + xx(or yy) are kinematically allowed for any given 
i, Dalitz-like plots will have 'stripes' extending from each of the endpoints of these 
decays to zero (or to m Xty if these are not approximately massless); these stripes will 
overlay the basic box/wedge structure outlined above. 

If the types of decay chains the Xi follow are known and in particular if one 
type dominates (e.g., two-body decays through one or more known intermediate 
states), the shape of the xx(yy) invariant mass spectrum can be predicted and this 
information used to compare densities of points in different regions of the Dalitz-like 
plot; this in turn allows one to measure ratios of cross section x BR for the different 
modes XiXj which are responsible for the various Dalitz shapes. The Dalitz-like plots 
can then provide information about dynamics in addition to kinematics (contained 
in the location of the endpoints). 

In any particular model there will be a set number of heavy particles expected 
to be produced at LHC energies; therefore the types of possible boxes and wedges is 
likewise set and the number of possible box-wedge combinations (with possible over- 
laying stripes) is fixed. Only some of these combinations are topologically distinct. 
For example, consider a sample of events where two XjXj-type production modes 
dominate. This will yield a Dalitz-like plot that looks like a 'box within a box' (ig- 
nore stripes for the moment). The topology alone would not indicate whether X 2 X 2 
and X3X3 are being produced or X 4 X 4 and X 5 X 5 are being produced. A collection of 
Dalitz patterns form a topological class if they can be transformed into each other by 
any amount of dilation; i.e., they can be deformed into each other without crossing 
any kinematical hard edges. Furthermore, a wedge of type ij is difficult to distin- 
guish from an ij- wedge and an w-box (Mjq < Mj.) combination and thus these two 
cases will be treated as topologically equivalent. If stripes are present the amount of 
degeneracy escalates (see subsequent application to SUSY). 

Application to Sparticle Decays 

The -R-parity conserving MSSM is next considered as a test-case for this technique. In 
the MSSM, there are four distinct neutralinos 3 , the lightest of which, Xi, is supposed 
to be stable and undetectable (e.g., in minimal supergravity-inspired models). These 
-inos are the physical eigenstates resulting from the two pairs of neutral electroweak 
(EW) gauginos and Higgsinos in the MSSM. Consider -ino pair production. More 
specifically, the modes of interest here are pair production of heavy -inos: 

pp -> intermediate(s) -> x°iX°j (hj = 2, 3, 4), (3) 

where both of the -inos subsequently decay (in the detector, as expected in all MSSM 
scenarios) leading to final states of the type described above. Thus neither -ino is 

3 In what follows, we will often refer to neutralinos collectively by the shorthand "-inos". 



Figure 1: Feynman diagrams for heavy = 2,3,4) neutralino pair production 
mechanisms: (a) 'direct' production via EW gauge boson; (b) Higgs- mediated production; 
and (c) production via cascade decays of gluinos (shown here) or via gluino/squark or 
squarks — to obtain these diagrams, make one or both of the squarks in (c) on-mass shell 
and remove the associated gluino(s) and the connected quark(s). 

allowed to be the stable Lightest Supersymmetric Particle (LSP) Xv However, in 
each event two LSP xV s are subsequently produced from the decays of the two initial 
heavy -inos, possibly through a chain of decays, along with the e + e~ and pairs 
we demand 4 . At the LHC, heavy -ino pair production occurs via virtual SM gauge 
bosons (termed "direct production"), via the decays of heavy Higgs bosons, or via 
cascade decays of coloured squarks and gluinos (see Figure HJ . The last of these, 
production via coloured intermediates, 

PP 99, M, QQ XiX°j = 2, 3, 4), (4) 

will be the focus of the current work 5 . Due to the strong coupling, (J1J) has the 
potential for yielding the largest number of signal events, if the intermediate gluinos 
and squarks are sufficiently light. Rates from EW direct production are typically 
too low for the technique described herein to be effectively utilised, except perhaps 
in certain minor regions of un-excluded parameter space or with allowance for ample 
time to gather more events Since the lightest MSSM Higgs boson (h°) can only 
yield LSP-containing -ino pairs, the Higgs-mediated production modes of interest 
involve the heavier MSSM Higgs bosons (H° and A ). Their masses need to be 
in the correct range to get sufficient Higgs boson production and yet have open 
decay modes to exclusively-heavy -ino pairs. This is certainly possible, as will be 
documented in another work [3J |H] . However, the E W production rate will lead to a 
smaller number of potential signal events than for (jlj. Thus, (@J should be the main 
source of -inos at the LHC if gluinos are light (~ 400 GeV). In the current work, 
inputs for gluinos and squarks will be set near the lower end of their allowed mass 
ranges while the input Higgs boson mass will be fixed fairly high up (~700 GeV) 6 . 

Aside from the larger possible signal rates with (J3| as compared to with the two 
EW production mechanisms, there are two other seminal distinctions between (J3J) and 

4 Other SM fermions aside from isolated electrons/muons may also be present in the final state. 
5 Here taken to include gx,qx production modes, which in fact only make minor contributions. 
6 These restrictions are reversed in a detailed look at the decays of heavy MSSM Higgs bosons 
into -inos in |S] . 



the other two that can strongly influence the analysis. First, as a side-product to 
producing -inos, the decaying gluinos/ squarks in (j3J) also typically lead to jet activity 
in the final state, whereas the other two production mechanisms may be hadronically 
quiet much or at least some of the time. Thus backgrounds to (@J) may be more severe 
/ less amenable to cuts. This could bring the signal rate after cuts down to the level 
of the other two processes. In fact, it will be shown in the simulations section to 
follow that the backgrounds are not so severe. Further, demanding the presence of 
jets is actually useful in reducing some backgrounds. 

The second point warranting attention is that in (0J) the two -inos are produced 
separately, whereas in the two EW processes there is an -ino-ino < -'- ) -S' vertex (where 
S is a SM gauge boson or an MSSM heavy Higgs boson). If the cascade decays 
were solely from gg and/or qq production, where here both squarks are of the same 
flavor and had the same SU (2) L quantum number (e.g., ul,u* l dRd* R , etc.), this would 
reduce the number of possible topologies that can result from (£Q) relative to the 
other production mechanisms (that is, considering XiXii X^x'j an d xlx^j production 
with i ^ j, knowing two of the three rates would determine the remaining one). 
However, gq production is very significant 7 (in fact, when all gq combinations are 
added together, their combined rate is larger than either the gg rate or the combined 
qq^ production rate [5]). If all the different squarks always decayed into gluinos, the 
afore-mentioned reduction in possible topologies would still occur. Actually, for the 
MSSM parameters herein considered, the different squarks decay into gluinos with 
BRs ranging from ~40% to ~95% (save for stops, which cannot decay into gluinos 
and top quarks in the cases examined), and the remaining times decay directly into 
charginos and neutralinos with differing BRs into the individual -inos, which would 
tend to restore the more general range of topologies if the different squark flavours 
contributed comparably. However, this is not the case — contributions from the ul 
squark are fairly dominant, and the BRs for this squark tend not to differ markedly 
from those of the gluino. So the reduction in topologies is partially true. How much 
this is so will be quantified later when specific points in the MSSM parameter space 
are discussed (see Tables |2] and EJ) • 

The possible Dalitz topologies from -ino decays to lepton pairs in (J2J) are built 
from 3 possible boxes and 3 possible wedges taken individually, and hence 63 basic 
combinations (Ef =1 C(6, i)) of boxes and wedges when considered all together, though 
of these many are topologically equivalent — only the 9 topologically distinct pat- 
terns shown in the left square of Figure 121 are possible. The patterns can be profitably 
categorised by the outer envelope exhibited (A, B, or C as shown in the right square 
of Figure |2J). Additional internal structure can then further sub-divide members of 
each envelope-type. To the extent that the reduction discussed in the preceding para- 
graph is applicable for (@J), the envelope- type then depends on the relative individual 
production times leptonic decay rates for X234 ( ca ^ these ^2,3,4) • If ?~4 is appreciable 
from both parent coloured sparticles, then box A is obtained (if r2 and/or are 

Further, squark-initiated processes are likely to contain extra jets which can increase the per- 
centage of these events that will pass a cut on the minimum number of jets in an event that will be 
imposed. 







Figure 2: At left, all possible distinct Dalitz patterns initiated by production of x^X^ pairs 
= 2, 3, 4) in the MSSM — 9 basic types exist which may be augmented by up to three 
stripes. At right, the region of the Dalitz-like plot can be broken up into the sections labeled 
a, (3, k. Envelope-types A, B, and C are possible. If the reduction from considering -inos 
produced in pairs to considering only individual -ino production rates is valid, only the 
patterns shown here are possible. The density of points inside each region of the Dalitz-like 
plot is shown to be uniform for simplicity only. 



also sizable, boxes and wedges inside of the box A envelope are also present); with 
r 2 > r 4 ^> r 3 , regions 5, rj, and k of the Dalitz plot (see Figure |2J) are down in popu- 
lation density by ~r 4 2 , and thus negligibly populated — resulting in wedge B; finally 
with r 2 ~ r 3 > r 4 regions rj and k again are negligibly populated yielding pattern C. 
Higgs-mediated -ino pair production would move beyond individual -ino production 
rates and probe the fundamental -ino-ino^-Higgs vertices, perhaps even more fertile 
subject-matter vis-a-vis application of the Dalitz-like technique [Hj (despite the lower 
maximal rates attainable). 

Most previous LHC studies [Zj of multi-lepton signals from gluino/squark cascade 
decays have concentrated on discovering evidence for SUSY, not upon extracting 
information about the sparticle spectrum from observed leptons' momenta. Ear- 
lier attempts to look at the spectra of invariant masses for lepton pairs resulting 
from gluino cascade decays are presented in [SJ El US UJ]- This work only exam- 
ined one-dimensional invariant mass spectra where the electrons and muons were not 
distinguished. Further, the work was restricted 8 (unlike the discovery search just 
mentioned) to the pair production of only the second lightest neutralino, X2%2- Sim- 
ilar restrictions are found in previous studies of Higgs-mediated -ino pair production 
[131 E|- In fact, in [Hj a Dalitz-like plot was presented, but the X2X2 _om y condition 
meant that only a box was possible. Thus the current work is novel for its inclusion 

8 52 also tried to obtain information on the -ino mass spectrum from a similar invariant mass 
reconstruction of tau-lepton pairs copiously produced in gluino decays at very high tan f3. Here 
mention is made of the heavier -inos in addition to 



of the heavier -inos, X3 an d X°ai the presence of which leads to a far richer variety of 
possible decay topologies for study via the Dalitz-like method. In addition, inclusion 
of the heaviest -ino states makes it more comfortable to construct sparticle spectra 
with slepton masses near or even below the heavier -inos. Such a sparticle mass 
hierarchy can greatly enhance the leptonic decay modes of the -inos - leading 
to far larger signal event rates. 

Decays of an -ino into a pair of same-flavour, oppositely-charged leptons plus the 
LSP may proceed through either two- or three-body processes with gauge boson or 
slepton intermediates; i.e., 

or ^-^i + ^r^-^+r + x? ' (6) 

where the two-body decays occur through an on-mass-shell Z°-boson or slepton and 
the three-body decays occur when the Z°-boson or slepton are off-mass-shell. The 
M(£ + £~) spectra from the different decay processes differ markedly. If the decay 
is via an on-shell Z°, then the lepton pair reconstructs the Z° and the spectrum 
is a sharp spike at Mz- If the decay is via an on-shell slepton, then the M(£ + £~) 
spectrum is basically triangular with a sharp rise in the number of event culminating 
at 



M(£+r)<m~ 1- pL L 



/ m~ \ 

1- — • 7 

V mj J 



Finally, if the decay is a three-body one via an off-shell Z° or slepton, then the 
M(£ + £~) spectrum is less sharply peaked toward the high end, but extends up to 

M{£ + £-) < m~o - m~o . (8) 

Regardless of whether the heavy states decay through two- or three-body de- 
cays, the distribution of dilepton invariant masses will be roughly triangular (i.e., 
more decays occur toward the endpoint). This has two immediate consequences: 1) 
hard kinematical edges in the Dalitz-like plot should be easy to identify since more 
of the event distribution is pushed up against the endpoint; 2) the distribution of 
points inside the boxes and wedges will not be uniform but can be fitted against 
an appropriate combination of triangular distributions: hence ratios of different -ino 
production cross sections contributing to an observed topology may be determined 
by comparing the number of points in different regions of the Dalitz-like plot. 

On the other hand, one or both of the initially-produced heavy -inos may not 
decay directly into the LSP plus leptons. Cascading decay chains including 

x2 - xl + tr , xl - x°2 + tr & x§ - xl + ?r (9) 

decays are also possible. If present, these will add the afore-mentioned stripes to the 
Dalitz-like plots. Note that up to three stripes are possible (these augment the 9 basic 
patterns shown in the left square of Figure |2J) • Figure El illustrates how the 9 basic 
patterns of Figure 121 may be altered by the presence of stripes to further enrich the 



Figure 3: Examples of how adding stripes affects the topologies shown in the left square 
of Figure El Some examples of one stripe added a) one stripe added to basic type #6, b) 
one stripe added to #2, c) one stripe added to #3, and d) two stripes added to #3 . 

number of possible Dalitz-like plot topologies. Note that among the topologies shown 
here, only in the case of Figure can the observed hard edges in the Dalitz-like plot 
be unambiguously linked to specific -ino -ino pairs including all the decaying -inos 
in the MSSM: the three boxes must correspond to X2X2, X3X3; and X4X4 modes, while 
the stripe must correspond to xt cascading through one of the other -inos (whether 
or not this is adequate to reconstruct all the mass differences in the complete -ino 
mass spectrum depends on the roles played by the sleptons). One can imagine quite 
elaborate decay chains, with Xi ~ y X3 ~~ * X2 ~ y Xi f° r instance. However, such 
elaborate chains are very unlikely to emerge from any reasonable or even allowed 
choice of MSSM input parameters. Further, each step in such elaborate decay chains 
either produces extra visible particles in the final state or one must pay the price of 
the BR to neutrino-containing states. The latter tends to make the contribution from 
such channels insignificant, while the former, in addition to also being suppressed by 
the additional BRs, may also be cut (or enhanced) if extra restrictions are placed on 
the final state composition in addition to demanding an e + e~ pair and a / u + /i~ pair. 
Another caveat is that decays with extra missing energy (carried off by neutrinos, 
for example) or missed particles can further smear the endpoint. 

To the discussion of leptonic -ino decays must be added the caveat that there 
are other possible xi decay modes where each lepton pair (e + e~ or does not 

emerge from three-body or multiple two-body decays of a single -ino — exactly two 
leptons of different flavour can be obtained from the same -ino (as noted in an earlier 
footnote). For example, consider the following decay chain that includes a chargino 
intermediate: 

x^^e+v + ^^e+ve'-p' + xl . (10) 



It is also possible for all four leptons to come from one of the initial -inos while the 
other -ino yields no leptons. This can occur, for example, if one -ino decays via 



(ii) 



while the other decays as 



wxx or \j 



QQXi 



(12) 



Again though such channels are at least somewhat suppressed by the additional 
required BRs. 

Finally, note that if the decays proceed via a chain of two-body decays including 
an on-mass-shell slepton, then the edge positions of the topological shapes will depend 
on the mass of the slepton involved, as seen 9 in ( jTj l. and this may give rise to an 
asymmetry (noted for one-dimensional endpoints in [5|) in the Dalitz-like plot if the 
sleptons are not degenerate in mass: boxes will become rectangles and wedges will 
no longer be symmetric under the exchange of axes. Whether or not such deviations 
from boxes and wedges are discernible depends on the slepton mass splittings. The 
present work will not address this issue (degenerate or nearly-degenerate selectron 
and smuon masses will be assumed). Evidence for direct slepton pair production may 
either be useful in determining what is going in -ino-ino pair production processes if 
sleptons have such low masses as to be produced with sufficient rates, or sleptons may 
be more massive and thus have low direct pair production rates so that the -ino-ino 
event topologies and rates may shed light on the else- wise inaccessible slepton sector 
of the model. 

In this study, sleptons will be kept fairly light so as to enhance the leptonic 
decay modes of the -inos [T5] . The experimental limits from LEP on the slepton 
masses are [16] m^^)^] > 99(91) [85] GeV and > 43.7GeV. To try to avoid 
producing leptons that are too soft, the charged sleptons (and the sneutrinos) are 
set sufficiently higher in mass than the LSP, which is the lightest neutralino, Xi> as 
mentioned earlier. Soft SUSY-breaking inputs are further simplified by assuming a 
flavour-diagonal slepton sector with = m~ lL and vanishing trilinear 'A-terms'. 
This effectively reduces the slepton sector to one soft SUSY-breaking input mass 
(identified with = m^ R = mi ) in the analyses that follow. These choices may 
be sub-optimal, especially since -ino decay modes to sneutrinos (which only depend 
on mi) tend to be 'spoiler' modes most often yielding only neutrinos in their decays 
and no charged leptons, and nothing prevents choosing a more complex set of inputs 
for which a topological analysis may yield even more information. Since -ino decays 
to tau-leptons are generally not anywhere near as beneficial as are -ino decays to 
electrons or muons, it would be even better if the stau inputs were significantly 
above those of the first two generations, thus the soft stau mass inputs are somewhat 
arbitrarily fixed to be 100 GeV above the degenerate soft mass input chosen for the 
first two generations. 

9 Note that in JJJ nig is the physical slepton mass (or masses, if more than one intermediate 
slepton is possible), not the soft slepton mass input, , defined in the following paragraph. The 
physical mass and the soft input may differ by several GeV or so. 



Table 1: Relevant sparticle masses (in GeV) for Points A, B and C. 





A 


B 


C 


Xi 


93.9 


186.4 


113.4 


xl 


136.7 


248.8 


175.8 


xl 


167.5 


257.8 


209.4 


xl 


236.9 


422.6 


295.2 


xt 


136.4 


238.3 


168.2 


xi 


238.2 


422.6 


295.1 


m~ 


136.60 


189.38 


135.47 




155.82 


204.74 


156.28 


m m 


155.75 


203.76 


154.67 


m eo 


156.71 


205.47 


157.24 




156.78 


206.44 


158.83 


m~ - m~ 


0.89 


0.73 


0.96 




1.03 


2.68 


4.16 



Since the soft mass inputs for selectrons and smuons are degenerate, the masses 
of the actual physical sleptons will also be nearly degenerate. With = m^ R and 
An = 0, the physical slepton masses are given by 



-Mf cos 2/3 ± 



m e R + m e 
m] R + \M 2 z cos 2/3. 



i (4M^ - 3M|) 2 cos 2 2/3 + m 2 /i 2 tan 2 /3 



(13) 
(14) 



The level at which the degeneracy is broken will be shown in some of the plots 
to follow; however, it remains too small to be quantitatively analysed. Thus the 
question of how non-degenerate selectron and smuon masses can affect the observed 
topologies will not be probed in this first realistic simulation. 



Detailed Study of Representative Points 

To illustrate this technique, three points in the MSSM parameter space with rep- 
resentative topologies were chosen for simulation utilizing the Monte Carlo package 
PH HERWIG 6.5. With common inputs of M A o = 700 GeV, m- g = 400 GeV, and 
rrig = 500 GeV (for all soft squark mass inputs), the three points are 

• Point A: tan/3 = 5, M 2 (M 1 ) = 200(100) GeV, /i = -150 GeV, m~ lR = 150 GeV. 

• Point B: tan/3 = 20, M 2 (M 1 ) = 400(200) GeV, /i = -250 GeV, m ia = 200 GeV. 

• Point C: tan/3 = 30, M 2 (M 1 ) = 250(125) GeV, /i = 200 GeV, m~ ijt = 150 GeV. 
Note that gaugino unification at a high (GUT) scale is assumed for the EW gaugi- 
nos, so Mi is not independent of M 2 (Mi ~ \M 2 at the EW scale). However, this 
restriction is relaxed for the gluino mass, which is taken as independent of the EW 



gaugino masses. The sparticle mass spectrum for these three points is shown in Ta- 
ble n Masses do not include radiative corrections, which are generally small. The 
shown physical slepton masses are also salient numbers to the Dalitz plot analyses as 
they enter into (JJJ) and also largely control the leptonic BRs of the -inos. Note that 
if left-right sfermion mixing — the m|/i 2 tan 2 (3 term in Eqn. (JT3Jl — is neglected, 
the mass splitting of the smuons becomes equal to that of the sleptons, and thus is 
evidently sometimes markedly under-estimated. Unfortunately, the physical slepton 
masses input (via IS ASUS Y 7.58 [TH]) into the HERWIG simulations do neglect this 
mixing. This will be seen in the Dalitz-like plots shown later. 

Table El gives the BRs for squarks and gluinos to decay into charginos and neu- 
tralinos, and for charginos and neutralinos to decay into final states with any number 
ofleptons. These BRs were calculated using IS A JET (IS ASUS Y) 7.58 HH] 10 - Naively, 
one might expect neutralinos (charginos) to only produce states with an even (odd) 
number of charged leptons. This is incorrect since combinations of quarks may also 
be produced in the decay chains, and said quark combinations can have a non-zero 
net charge. Note the significant BRs for xf — > 2£ in Table EJ As Figure 1(c) clearly 
shows, quarks are expected even before neutralino/chargino decays are considered - 
demanding hadronically-quiet events is not an option in this case (but may be with 
the other production modes), in fact just the opposite is most effective: a minimum 
jet requirement will in fact be employed in the analysis to follow. The neutralino 
and chargino BRs to nC. (n = 0,1,...) final states given in Table El include neither 
r- leptons nor is from r-decays. Also, no demands are made on the leptonic px or \r)\ 
values. 

Combining the leptonic BRs of the assorted neutralinos and charginos with their 
production rates from decays of gluinos and the different squarks yields Table El 
The neutralino, chargino, or neutralino/chargino pair listed represents the first EW 
sparticles produced in decays of the colored sparticles. The EW sparticles can then 
themselves decay into other EW sparticles. 'Mixed' production modes (gXii QXii 
9xti QXi) are a l so included. These mixed modes account for only about (2.7%, 1.0%, 
1.4%) of the events for MSSM Point (A, B, C). Since there is less jet activity with the 
mixed modes, they are more likely to fail the minimum jet requirement. HERWIG 
lacks the facilities for giving the cross-sections for each separate gq, qcf'\ gx and qx 
process, so these were calculated using ISA JET 7.67^1 with CTEQ5jTH] parton dis- 
tributions. Some tinkering with the HERWIG code was however able to yield values 
for a(J2 M) an d cr(J2 qq^) as well as cr(gg). ISAJET+CTEQ5 cross-sections were vir- 
tually always found to be lower than those from HERWIG+CTEQ6. The <t(YJ qq^)s 
agreed to ~5% at the three MSSM parameter points, while the ISAJET+CTEQ5 
o-(E<?#)s HEff?» were lower by roughly 10-20% (5-10%). Given that HERWIG 
and ISAJET differ in the scales adopted for the parton distribution functions (which 
are also different here) and for the evolution of coupling constants, the differences 
seen in these cross-sections are in fact quite modest. Thus using ISAJET rather 
than HERWIG values should not markedly effect the estimates obtained. For an 

10 This is the version incorporated into HERWIG 6.5; however, results sometimes differ significantly 
from those obtained with later versions of ISASUSY. 



Table 2: Chargino and neutralino production and decay properties for Points A, B and C 
obtained from ISAJET(ISASUSY) 7.58: gluino & squark BRs to neutralino/chargino and 
the inclusive BR for obtaining any number of leptons (£ = e, /i) from a neutralino/chargino 
(includes neither r- leptons nor leptons from r-decays) . Squark decays include decays via a 
gluino; here X denotes some number of quarks and gluons. Reverse the signs given for the 
charginos for anti-squark decays. "— " means no BR at tree level. 



BRs for 


Point A 


Point B 


Point C 


g ->_xi g/qq 

u L /d L /u R /d R -> Xi x 


.23 

.11/.13/.34/.26 


.61 

.49/.54/.64/.62 


.34 

.18/.22/.43/.36 


9 g/qq 

u L /d L /u R /d R -> xl x 


.19 

.15/.13/.19/.19 


.18 

.18/.16/.18/.18 


.25 

.21/.19/.24/.25 


Xl - M/U/2Z 


.76/~10~ 7 /.25 


.001/-/.999 


.305/-/.695 


9 ^xl g/qq 

u L /d L /u R /d R -> xl x 


.038 

.02/.03/.03/.04 


.11 

.09/.09/.09/.10 


.10 

.05/.06/.08/.09 


Xl - M/U/2t 
/M/U 


.93/.004/.07 
/~10" 9 /.002 


.76/.00045/.24 

/-/- 


.80/.005/.19 
/-/3 x 10~ 5 


g ^xl g/qq 

u L /d L /u R /d R -> 


.11 

.15/.17/.09/.10 




.03/.03/~io- 5 /-io- 5 


.024 

.09/.10/.02/.02 


x2 - 0£/l£/2£ 
/M/U 
jUjU 


.63/.094/.26 
/~10- 8 /.019 
/-/- 


.46/.16/.38 
/.0004 /.0013 

/^io- 11 /~io- 7 


.47/.20/.33 
/.00008 /.0004 

/- /~io- 7 


g -^Xi qq' _ 

u L /d L /u R /d R -> x^X 
u L /d L /u R /d R -> xT^ 


.28 

.25/.07/.11/.13 
.07/.12/.11/.13 


.11 

.12/.05/.05/.05 
.04/.07/.05/.05 


.26 

.25/.07/.10/.12 
.07/.16/.10/.12 


x? - 0£/H 


.61/.39 


.0023/.9977 


.0016/.9984 


u L /d L /u R /d R -> x^X 
u L /d L /u R /d R -> X2^ 


.16 

.21/.04/.06/.07 
.04/.30/.06/.07 




.05/-/-/- 
-/.07/-/- 


.033 

.14/.010/.013/.016 
.009/.20/.013/.016 


/2£/3£ 
/4£/5£ 


.11/.83 
/.020 /.038 
/^10- 6 /~10" 7 


.135/.756 
/.084/.025 

/~io- 9 /~io- 6 


.060/.853 
/.065/.021 

/^10- 7 /~lQ- 7 



integrated luminosity of 30 fb 1 (equivalent to two or three years of low-luminosity 
performance at the LHC) ISAJET+CTEQ5 predicts approximately 60,000, 200,000 
and 197,000 4£ events before any cuts are applied for MSSM Points A, B and C, 
respectively. By contrast there are only 744, 471 and 750 4£ events from 'direct' 
production of charginos/neutralinos at the three points in parameter space (1.2%, 
0.2% and 0.4% of the coloured-sparticle cascade rates). 

The percentages given in parentheses in Table E3 are when only production via 
gluinos is considered. Here the reduction in the possible topologies mentioned ear- 



Table 3: Percentage contributions to 4£ events from the various neutralino/chargino 
pair production modes for MSSM Points A, B and C. Numbers in parentheses only 
consider gluino pair production. 



Point A 


Point B 


Point C 


,rO,rO 

X2X4 


OT A C 0/ 

27.45/0 


/OK 00/ \ 

(25.2/o) 


X2X2 


to oy 


(76.37o) 


X2X2 


66.6/0 


/TO TO/ \ 

(72.7/o) 


~0 ~0 

X2X2 


14.5% 


/1 r T07 \ 

(16.7%) 


~o ~o 
X2X3 


20.3% 


/oo 1 0/ \ 

(22.1%) 


~o ~o 
X2XI 


1 00/ 
13.3% 


(15.9%) 


X4X4 


11/1 K®f 
11.40/0 


{0.0/0) 


X2X4 


Q 9O/ 


(n r\<Y\ 
(U.U/o) 


X2X4 


1Z.D /o 


/<? C0/\ 

(0.0/0) 


xlxl 


7.8% 


(7.6%) 


xlxl 


1.4% 


(1.6%) 


xixi 


2.7% 


(1.8%) 


X2X2 


7.4% 


(6.25%) 


xixi 


0.9% 


(0.0%) 


xlxl 


1.3% 


(0.7%) 


xixi 


6.6% 


(8.7%) 


xlxl 


0.5% 


(0.0%) 


X1X2 


0.74% 


(0.9%) 


xixi 


6.5% 


(6.25%) 


xtxi 


0.1% 


(0.0%) 


xlxl 


0.68% 


(0.9%) 


xixi 


4.4% 


(5.15%) 


xtx2 


0.1% 


(0.0%) 


xixi 


0.67% 


(0.9%) 


xixi 


4.3% 


(5.15%) 


xtxi 


0.08% 


(0.0%) 


xlxl 


0.5% 


(0.15%) 


xixi 


2.9% 


(2.9%) 








xixi 


0.3% 


(0.2%) 


xlxl 


2.6% 


(2.4%) 








xixi 


0.3% 


(0.1%) 


xixi 


1.9% 


(2.3%) 








X2X2 


0.24% 


(0.05%) 


xlxl 


1.8% 


(2.1%) 








xixi 


0.15% 


(0.05%) 


xMxl 


0.2% 


(0.3%) 














xixi 


0.1% 


(0.2%) 














xixi 


0.1% 


(0.2%) 















lier 11 applies. For instance, look at the xlxl^ X3X3 an d xlxl fractions (in parentheses) 
for Points B and C — or the xlxl, xlxl an d xlxl fractions for Point A — labeling 
these as rf, r| and r^-, respectively (they are proportional to the production cross- 
section times the 4 lepton BR for the given -ino pair), we find that r\j = 2r\rj. 
Explicitly (ignoring the insignificant 'mixed' production channels), 

BR(2 coloured sparticles — > neutralinoj + neutralino.,) 

* BR(neutralinOj + neutralinoj — > 41eptons), (15) 

BR(gluino or squark — > neutralinoj) 

* BR(neutralinoj — > designated number of leptons). (16) 

Under the assumptions that 

• 1. Each -ino came from a gluino; 

• 2. Each -ino produced two leptons 
then Tij factorises as 

r\j = BR(2 gluinos — > neutralinoj + neutralino,,) 

* BR(neutralinOj + neutralino^ — > 4 leptons 

11 The so-called 'second point warranting attention' in the section entitled "Application to Sparticle 
Decays." 



while 



(presumably of the same flavour) ; 



2 * BR(gluino — > neutralinoj) * BR(neutralhiOj — > 2 leptons) 
* BR(gluino — ► neutralinoj) * BR(neutralinOj — > 2 leptons) 
2nrj . 



(17) 



Checking the BRs for the various -ino pairs from when one -ino produces m leptons 
and the other produces n leptons, where m + n = 4, shows that only the m = 
n = 2 case contributes significantly. The difference between the percentage when all 
production modes are included and the percentage in parentheses thus quantifies the 
deviation due to the squark production modes. That these two values generally do 
not differ by too much indicates that the relationships among the -ino pair production 
rates expected from gluino-only production do to a significant extent remain intact 
when squarks are included. There is a caveat to this though: here only inclusive 4£ 
events are tabulated with no cuts; squark events may contain more jets and thus a 
higher percentage of them may pass a minimum jet number requirement. 

Consider a numerical example (to be compared later to results extracted from 
simulations via the Dalitz-like technique): for Point C, Table El gives 



(or r 23 : r 24 : r 34 = .159 : .0660 : .00721 = 22.0 : 9.16 : 1 if only gluino pair 
production is considered). Assuming the formula = 2r^j holds, it follows that 
r 2 : r 3 : r 4 = 19.00 : 1.04 : 1 (r 2 : r 3 : r 4 = 22.00 : 2.404 : 1 with only gluino pair pro- 
duction 12 — this result matches the values obtainable from Table |2J). Alternatively, 
the identical -ino pair values can be used, assuming r# = r\. For Point C, Table El 
then gives 



(or r 22 : r 33 : r 44 = .727 : .009 : .0015 = 485.3 : 5.78 : 1 if only gluino pair production 
is considered). These values yield r 2 : r 3 : r 4 = 11.47 : 1.16 : 1 (r 2 : r 3 : r 4 = 22.03 : 
2.404 : 1 with only gluino pair production). Note that the results considering only 
gluino pair production agree, while the full results do not. Thus disagreement in 
such calculations indicates significant contributions from squark production. 

Roughly a third of the 4£ events for Point A come from production modes includ- 
ing charginos. However, a substantial fraction of these events will not have leptons in 
same-flavour, opposite-sign pairs. So their effect on this analysis will be diminished 13 . 
This does expose a minor weakness of the framework developed herein which is built 
only for the neutralinos. The chargino production contributions for Points B and C 
are much smaller (~1% and ~5%, respectively). 

Table El is only expected to serve as a guideline against which simulation results 
may be examined. While this will prove useful in confirming the interpretations of 

12 The ratio of <?g-production to gq- and qif 1 ' -production is crucial here, and this ratio is larger 
for HERWIG+CTEQ6 than for ISAJET+CTEQ5. So the latter would yield larger deviations from 
Tij = 2riTj. For such deviations to be taken as evidence for squark-initiated processes, it needs 
to be shown that the measured deviations exceed the uncertainties due to structure functions and 
simulator cross-section estimates. 

13 In fact, same-flavour, like-sign 4£ events (i.e., e ± e ± /^ =F /j, =F events) could be used to estimate the 
chargino contribution and then remove it. This is seen, though with a different rationale, in |9lllU|. 



r 23 : r 24 : r 34 = .133 : .125 : .013 = 10.2 : 9.6 : 1 



(18) 



r 22 : r 33 : r 44 = .666 : .0068 : .005 = 131.5 : 1.3 : 1 



(19) 



features on the Dalitz-like plots, it should be emphasised that this is information 
that the real experiments will not be able to access; i.e., with a simulation, we are 
of course able to choose what point in the parameter space to simulate. 

Numerical Simulations 

Events for pp —> gg, gq, qq, gx, qx were generated at the specific points in the MSSM 
parameter space just discussed using HERWIG 6.5 coupled with a detector simulation 
which assumes a typical LHC experiment, as provided by private programs checked 
against results in the literature. The CTEQ6M JH] set of structure functions was 
used in conjunction with HERWIG 6.5 to determine the cross-sections. 

In event selection and subsequent cuts, stress is put on keeping the cuts reason- 
ably general so that they will hopefully be applicable across a large swath of the 
allowable parameter space. These cuts can quite probably be further honed once the 
first evidence(hints) of possible MSSM events is discerned, and the rather minimal 
selection criteria used here certainly do not represent the optimal choice for the few 
points examined in detail in this work. The actual criteria used are as follows: 

• 4£ events: events are selected which have exactly four isolated leptons 

(£ = e ± or / u ± ) with \rf\ < 2.4 and > 7, 4GeV for e ± 1 n ± 1 respectively. 
The isolation criterion demands there be no tracks (of charged particles) with 
Pt > 1.5 GeV in a cone of r = 0.3 radians around a specific lepton, and also 
that the energy deposited in the electromagnetic calorimeter be less than 3 GeV 
for 0.05 radians < r < 0.3 radians. 

• The four leptons must consist of exactly one e + e~ pair and one n + fi~ pair. 

• A cut on missing transverse energy: events must have Ifi? > 20 GeV. 

• Three or more jets must be present. Jets are defined by a cone algorithm with 
r = 0.4 and must have \rf\ < 2.4 and E 3 T > 20 GeV. 

In addition to requiring the advertised lepton make-up of the final state demanded by 
the Dalitz-like technique, we apply the well-known canonical missing energy cut to 
select for SUSY events. Further, we require some minimum number of jets as noted. 
Since production mechanism (J3J) typically generates considerable jet activity, while 
the H, A 'background' does not, such a cut is found to be particularly helpful. There 
are no cuts on the momenta properties of the leptons (aside from demanding that they 
be hard and central enough). This is consistent both with (a) the wish to paint all the 
leptons from the signal events onto our Dalitz-esque canvass to show the richness of 
the possible topologies, and with (b) the desire not to narrow the applicability down 
to only a minor portion of the MSSM parameter space currently experimentally 
viable that can satisfy such additional restrictions. The lepton isolation criterion is 
essential though to virtually eliminate enormous QCD backgrounds from events with 
leptonically-decaying 6-quarks (such as from ti production). Note also that gluino 
cascade decays are often rich in 6-quarks (particularly for higher values of tan /3); thus 



Table 4: Number of events after the successive cuts denned in the text for MSSM Parameter 
Points A, B, and C (for 30 ftr 1 ). 





Process 


A£ events 


e~*~e -jj^u pairs 


rp miss 


AT ^> Q 
J V jets ^ o 




Z°Z° 


365 


175 


11 


o 




Z° + jet 














SM processes 


t H + , t H~ 


1 


1 


1 


o 


(common) 


ti 
















ttz° 


47 


7 


6 


2 




tth° 


4 


1 


1 







total SM bkg. 


417 


184 


19 


2 




y y (direct) 


26 


6 


6 





Point A 


I,v 


50 


23 


21 


1 




A°,H° 


29 


5 


5 







g, q signal 


7628 


2350 


2292 


2110 




XX (direct) 


95 


36 


34 


2 


Point B 


i,v 


5 













A°,H° 


85 


26 


26 


4 




g, q signal 


15652 


7114 


6883 


5979 




XX (direct) 


78 


20 


18 


1 


Point C 


£,u 


12 


3 


3 







A°,H° 


292 


114 


109 


5 




g, q signal 


19897 


8595 


8357 


7615 



addition of a 6-tagging requirement might further enhance the signal/background 
ratio — but probably at the loss of some significant fraction of the signal events. For 
present purposes such a cut was found unnecessary. That, as will be shown next, the 
signal stands out over the backgrounds with so few cuts attests to the robustness of 
this signature and to the potential to obtain Dalitz-like plots using realistic simulated 
data that reflect the theoretical expectations discussed in the previous sections. 

Both g,q signal events and SM & MSSM background events were simulated at 
each point again assuming an integrated luminosity of 30 ftT 1 to see how the expected 
features show up for a realistic sample size. The number of events passing each of the 
cuts above for points A, B and C are listed in Table EJ which clearly illustrates how 
effective even this limited set of cuts is at eliminating the SM backgrounds 14 . And, 
at least for the particular cases studied here in detail, a sufficient percentage of signal 
events pass the cuts, while the number of surviving events from SUSY "background" 
processes — primarily pp Xi' Xi^i Xi^Xj' an d PP — > A , H° — is negligible. 

14 SM background events from Z°Z° production would be concentrated around Mz were they 
not eliminated by the three jet minimum requirement. When the other production mechanisms are 
considered in later works, such a requirement will probably neither be possible nor desirable. Then, 
though the relative number of SM background events passing cuts may yet be quite low relative to 
the total number of signal events in the plot, they can still lead to uncertainty in precisely locating 
edges (particularly indistinct ones) that happen to be in the close vicinity of Mz- 



Note that the '4£' rates given in the first column of Table 0] are roughly an 
order of magnitude smaller than the inclusive 4£ rates predicted in the previous 
section from the ISASUSY inputs. Further investigation indicated that somewhat 
less than half of the events were lost when the lowest Et lepton failed the imposed 
minimum Et cut. Other factors in the event (such as heavy quark decays) could also 
have yielded extra leptons, so that 4-lepton events became 5-lepton events; however, 
the number of n > 4 lepton events was checked to be quite small. Some events 
certainly had leptons too close to the beam pipe, but, again, this is not expected to 
be a major factor. We are thus led to conclude that the majority of the 4£ events 
were removed due to the isolation requirements. The fact that, as we shall see, 
the simulation results, qualitatively at least, track the values given in Table |3] fairly 
well is consistent with this hypothesis (if the main factor had been the minimum 
Et cut, for instance, XiX*j events, where i and/or j is 2, might have been highly 
preferentially eliminated). Nonetheless, the large fraction of events removed and 
the subsequent cuts applied caution against expecting a high degree of quantitative 
agreement between the simulation results and those of Table El (as already noted). 

Figure 0] presents the Dalitz-like plot for MSSM parameter point A. A 'wedge 
inside of a box' topological structure is apparent (as per pattern A in the right-side 
square of Figure EJ), a clear indication that two pairs of -inos, xlx^j (i < j) and 
are being produced at significant rates. The (j3J) production mechanism then demands 
that a XiXi box also be present, the position of which overlaps with that of the low 
M(e + e~), low M{n + ji~) corner of the wedge (as noted earlier, adding such a box is 
not viewed as being topologically distinct). A hard kinematical edge [i.e., the line 
in the plot across which the density of points changes very rapidly) at ~40-45 GeV 
is very apparent. The outer box seems to end at ~140 GeV though there are a small 
number of straggling points beyond this mostly at high M(e + e~), low M(/i + /i~) and 
at low M(e + e~), high M(/z + /i~). Also discernible inside the wedge are somewhat 
indistinct drops in population densities along both axes at ~85 GeV. 

The shaded bands and dashed lines included in the plot show the expected lo- 
cations of hard edges based on the -ino and slepton mass spectrum obtained from 
ISASUSY for Point A. The ~40-45 GeV hard edge corresponds to the 42.8 GeV 
XtXi mass difference. Here X2 is decaying through an off-shell Z° or slepton, with 
BR(x2 — > Xi^ + ^~) = 0.245 (very unlike the leptonic BR for the Z°) indicating that 
the off-shell sleptons are playing large roles. The other -inos decay mainly through 
on-mass-shell sleptons. The outer edge at ~140 GeV agrees with the endpoint for 
the two-body decay chain xt ~ * ^ — * Xi + ^ + ^~ (though this is the actual decay 
channel for this sparticle spectrum, in fact the two-body decay and three-body decay 
endpoints differ by less than 1 GeV in this case). So the outer box is from XaXa pro- 
duction and the wedge is from X2X4 (including an inner box from X2X2 production). 

The population changes at ~85 GeV inside the wedge might be interpreted as 
evidence for significant X2X3 production, or as a 'stripe'. In fact, they are due to the 
latter, and are associated with the decay chain xt — > H>& — ► X2 + ^ + ^~ which happens 
22.8% of the time. The xt decay chain mentioned in the last paragraph occurs 65.8% 
of the time, and the remaining 10.6% of the xt decays are through xt- Note that at 
least some a priori knowledge of the -ino mass spectrum and decay modes is required 




Figure 4: M(e + e~) vs. M(^i + fi~) Dalitz-like plot for MSSM Point A assuming an integrated 
luminosity of 30 fb -1 . Signal events from gluino/squark production are denoted by dots, 
SM background events by crosses and events from other SUSY processes by open circles. 
Shaded bands (dashed lines) indicate kinematical endpoints expected from two-body decays 
(three-body decays) based on ISASUSY. 



to designate this feature a stripe, showing that such Dalitz-like plots do not always 
uniquely identify the underlying -ino production/decay modes. Note also that the 
position of this feature is given by (J7J), with m~o replaced by m~o, which in this case 
is quite different from m~o — m~o = 100.2 GeV. Thus care must be taken before 
assuming that features in invariant mass plots correspond to -ino mass differences. 

The designations in the last two paragraphs agree well with the percentages given 
in Table 01 including the 'stripe' assignment above as well as the absence of a 
associated box or wedges in Figure 0] The events lying outside the outer box in the 
Dalitz-like plot are due at least in part to production modes including charginos. 
This was confirmed in the HERWIG simulation by checking the identities of the 
parent particles of the leptons in these outlying events. In addition, a sampling of 
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Figure 5: M(e + e~) vs. M(/i + /i~) Dalitz-like plot for MSSM Point B assuming an integrated 
luminosity of 30 fb -1 . Signal events from gluino/squark production are denoted by dots, 
SM background events by crosses and events from other SUSY processes by open circles. 
Shaded bands indicate kinematical endpoints expected based on ISASUSY. The insert in 
the upper right corner zooms in on a critical region, and the assumed integrated luminosity 
is raised to 60 fb -1 . 

such events were also found to have leptons from top-quark decays or lost leptons 
(i.e., they were 5 lepton events with one of the leptons being too soft to pass the 
minimum Et cut or too close to the beam axis). 

Figure 03 for Point B displays a somewhat sparsely-populated wedge envelope 
matching pattern B in the right-side square of Figure 121 The interior edges for the 
wedge are at ~60-65 GeV, and event points taper off around 140-180 GeV. Inside 
this wedge is a much more densely-populated box with edges at roughly 60 GeV. A 
second very short-legged wedge structure is also indicated, with edges at 60 GeV and 
65 GeV. More events from a longer run time would help to clarify the structure in the 
crucial (M(e + e"), M(/i+ /U -)) = (60-65 GeV, 60-65 GeV) region of the Dalitz-like plot 



(see insert in Figure EJ). The plot bespeaks of dominant XiXi production with weaker 
contributions from XiX'j an d xlxl < j < k) (the latter yielding the outer wedge 
envelope and the former the short-legged wedge). Note in the MSSM framework i, j 
and k must be 2, 3 and 4. The short, stubby wedge tells us that two of the heavier 
-inos, presumably x\ and x% are quite close in mass. This is in very good agreement 
with the predictions from Table El a densely-populated X2X2 box and a short, stubby 
X2X3 wedge (the X3X3 box is too sparsely populated to be recognised — again the 
(M (e+e - ), M{n + li-)) = (60-65 GeV, 60-65 GeV) region of the Dalitz-like plot is seen 
to be crucial, with more statistics desirable to clarify the situation. Also, for this 
point in MSSM parameter space, m~o is in fact rather close to Tn~o. 

Shaded bands in the plot again show the expected locations of hard edges based on 
the -ino and slepton mass spectrum obtained from ISASUSY. Though the 62.4 GeV 
XtXi mass difference from ISASUSY roughly fits the position of the box edges, 
ISASUSY also reveals that the X2 decays nearly always through an on-shell slepton, 
BR(%2 ~~ *■ M — > Xi + ^ + ^~) = 0.999, with the lighter (predominantly right) and 
heavier (predominantly left) slepton mass eigenstates contributing about equally. 
Significantly, the spoiler decay modes to sneutrinos only have a suprisingly low BR 
of only ~10 -3 . Applying (J7|) using only the physical selectron masses (to match the 
HERWIG inputs) from Tabled predicts edges at 58.5 GeV and 59.0 GeV, confirming 
that the inner box is from X2X2 production. Again, ^3 almost always decays via 
on-shell sleptons to Xi, but now 71.6% of the decays are into sneutrino spoiler modes 
yielding no charged leptons. Application of (J7J) now predicts endpoints at 64.8 GeV 
and 65.5 GeV, about 6 GeV less than m~o — m~ . 

Note that gluino decays to xl 01 xf are kinematically impossible and events 
including a gluino decay to xt cannot generate 41 events. Yet Table El says that 
3.0% of the 4£ events are from X2X4 production and an outer X2X4 wedge is clearly 
visible in the Dalitz-like plot. This outer wedge must be due solely to production of 
the heavier squarks which are heavy enough to allow decays to Xa (the very small 
contribution from X4X4 is insufficient to generate the apparently-missing outer box). 
The more massive xl decays 59.4% of the time into sleptons (26.6% of the time into 
charged sleptons and 32.8% of the time into sneutrinos). The predicted endpoints for 
the charged slepton decays from (0) are 153.0 GeV and 155.4 GeV, basically giving 
the outer ends of the wedge envelope over 80 GeV below the XatXi mass difference. 
Again, some if not all of the events lying outside these bounds come from processes 
involving charginos. 

Lastly, Figure |H1 shows the Dalitz-like plot obtained for MSSM parameter point 
C . A wedge extending out to ~175-180 GeV is readily seen. The inner edges of this 
wedge are at ~55 GeV. Inside of this wedge is a shorter wedge with the same inner 
edges terminating at ~95 GeV, and in the corner a densely-populated box. Structures 
outside this wedge are more difficult to discern with this number of events: a box 
with edges at 95 GeV is somewhat clear while a wedge with inner edges at ~95 GeV 
extending out to ends at ~175-180 GeV may be barely discernible. Thus this plot 
could be classified as either pattern B or pattern C according to the nomenclature 
introduced in the right-side square of Figure 121 The features exhibited suggest fairly 
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Figure 6: M (e + e~) vs. M(/i + /x~) Dalitz-like plot for MSSM Point C assuming an integrated 
luminosity of 30 fb _1 . Signal events from gluino/squark production are denoted by dots, SM 
background events by crosses and events from other SUSY processes by open circles. Shaded 
bands indicate kinematical endpoints expected based on ISASUSY. The band marking the 
endpoint for X3 ~~ > i — > decays is quite narrow, and is drawn as a dashed line in 

this plot. 

dominant XiXi production, but with significant contributions from XaX^ and xlxl 
(i < j < k), and with lesser but still detectable contributions from XjXj an d X*jXk- 
(Again, in the MSSM, i, j and k must be 2, 3 and 4.) 

ISASUSY numbers for the -ino and slepton mass spectrum again yield the shaded 
bands and dashed lines showing the expected locations of hard edges. X2 virtually 
always decays via on-shell sleptons to Xi, 69.5% of the time through charged sleptons 
and 30.5% of the time through sneutrinos (contrast this with X2 decays at MSSM 
parameter point B). Again applying (j2J) using only the physical selectron masses 
(to match the HERWIG inputs) from Table [T] leads to predicted edges at 54.4 GeV 
and 55.3 GeV, corroborating that the inner box is from X2X2 production. Note that 



use of more correct physical smuon masses incorporating left-right sfermion mixing 
would significantly widen the horizontal shaded bands in Figure El 

Comparing the wedges in Figure E] with the one in Figure El we conclude that in 
this case the -ino masses are not so close together, decays via on-shell sleptons 
77.6% of the time (17.7% via charged sleptons and 59.9% via sneutrinos), the rest 
of the time decaying via an on-shell Z° (21.6%) or a xt (0.7%) or a X2 (< 0.1%). 
The charged-slepton-mediated x% decays should have endpoints at 95.8 GeV and 
95.9 GeV (in this case the two-body endpoint is nearly equal to the X3-X1 mass 
difference). The X3 decays via Z° lead to a band at 91.2 GeV also faintly visible in 
Figure El As with MSSM parameter point B, the sneutrino spoiler modes are much 
stronger (considerably stronger) in X3 (xi) decays than for X2 decays, suppressing 
contributions from the former to the Dalitz-like plot relative to the latter. Decays of 
X4 via on-shell charged sleptons (which occurs 26.9% of the time, compared to 43.7% 
of the decays being via sneutrinos) will result in edges at 172.3 GeV and 173.0 GeV 
(10 GeV or so below m~o — m~o). Xa also decays 15 22.3% of the time into XiW T which 
can yield aberrant events not anticipated in the neutralinos-only framework followed 
here. The variation in the widths of the shaded bands due to decays occurring 
through the two different same-flavour sleptons, which are 4.06, 0.44, and 3.24 GeV, 
can readily be understood from the variation of (|7|): 



In this case the x% ~^ Xi an d Xa ~~ * Xi bands have similar widths, as this is inversely 
proportional to the endpoint yet partially compensated for by the factor in paren- 
theses for much heavier Xa'i this factor is however very small for intermediate-mass 

where m~o 2 m~o 2 /m/ « 1, hence the relatively thin xt —> Xi band. 

Summarizing, the predicted endpoints from charged-slepton-mediated decays of 
X3 and xt affirm that the shorter (longer) wedge is from X2X3 (X2X4) production. The 
more faintly discernible box with edges at ~95 GeV is attributed to X3X3 production, 
the even more faint wedge of which this box is the corner is from X3X4 production, 
and the few events in the upper-left corner of the plot are presumably from XaXa 
production. The relative percentages of 4£ events given in Table El agrees fairly well 
with the densities of points seen in the associated features in Figure El 

It is interesting to see how effectively the relative -ino pair contributions can be 
extracted from the Dalitz-like plot. Assuming some knowledge of the dilepton invari- 
ant mass distribution, an estimate of the ratio of the different -ino pair production 
rates (stemming from the gluino/squark BRs to the different -inos) is obtainable 
from counting the total number of points in each region of the Dalitz plot and then 
taking the ratio. Approximating the distributions as being exactly triangular (see 
P), and taking the endpoint locations noted in the preceding paragraphs as 55 GeV 

15 There are also a smattering of other Xa decay modes: to staus 1.4% of the time, to v T 3.8%, 
-> xlixl) + h° 1.2%(0.1%), and -> xUxf) + Z° 0.2%(0.1%). These could only contribute a very 
small fraction of the events. 




(20) 



, 96 GeV and 173 GeV, the following rate comparison can be extracted 



r 22 : r 23 : r 24 : r 33 : r 34 : r 44 = 431 : 118 : 59 : 15.5 : 9.4 : 1 , (21) 

where is the rate from xl Xj production, or, considering just the the three wedges, 
'■ ^ 2 4 : ^34 = 12.3 : 6.3 : 1 . Compare these values to the results obtained earlier 
from Table El 

r 22 : r 33 : r 44 = 131.5 : 1.3 : 1 and r 23 : r 24 : r 34 = 10.2 : 9.6 : 1 . (22) 

Little more than crude agreement is discernible; bear in mind though that, as noted 
earlier, discrepancies may be reasonably expected in comparing all-inclusive 4£ rates 
with e + e~yU + yU~ rates after cuts. The assumption of strictly triangular population 
profiles is also certainly somewhat inaccurate. And, at least with modest statistics 
(i.e., with only results from the first year or two of running for the LHC), there will 
be significant imprecision in pinpointing the locations of the endpoints (the main 
source of uncertainty in the calculation if the triangular distribution assumption is 
viable). One factor that is not an important concern at this MSSM parameter point 
is contamination from chargino-related events; but, said contamination could skew 
such a calculation at other MSSM points (as, for instance, Point A). 

If = 2r^j is now assumed to be valid (note though the afore-mentioned serious 
caveats to this assumption), the relative individual production times leptonic BRs 
are obtainable: r 2 : r 3 : r 4 = 12.3 : 1.6 : 1 (using r 23 , r 34 and r 22 or r 24 as inputs) 
or r 2 : r 3 : r 4 = 12.6 : 4.7 : 1 (using r 23 , r 34 and r 44 as inputs). The extent 
to which these two results disagree could (as before for the inclusive A£ results) be 
interpreted as implying significant contributions from squark-production (though the 
inaccuracy of the triangular distribution assumption may also be a factor). Recall 
inclusive 4£ results were r 2 : r 3 : r 4 = 19.00 : 1.04 : 1 (using r 23 , r 34 and r 24 
as inputs) and r 2 : r 3 : r 4 = 11.47 : 1.16 : 1 (using r 22 , r 33 and r 44 as inputs). 
Apparently, dynamical information from the densities of events in the Dalitz-like 
plot's various geometrical components may be more difficult to extract than the 
kinematical information contained in the location of the hard edges. However, more 
sophisticated statistical analyses may be expected to yield better results. 

Next contrast the information apparent in the Dalitz-like plots with that readily 
obtainable from the more traditional one-dimensional projections shown in Figure 
Notice how similar the results for Points B and C appear in Figure while Figure 
Eland Figure IH1 are quite different. Note also in this case the sharp drops observed 
would only be sufficient to identify which -inos are being produced, not which -ino 
pairs are being produced. 

Concluding Remarks 

Production of pairs of new heavy particle states XiXj at hadron colliders has been 
studied emphasising the simple topological forms expected in certain two-dimensional 

16 Calculational details are relegated to an appendix. 
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Figure 7: One-dimensional projection of Figures for MSSM Points A, B and C assuming 
an integrated luminosity of 30 fb -1 . Information about the decay topology is lost in such 
projections. 



Dalitz-like plots. It is assumed the heavy particles decay into pairs of SM particles 
(with the possible addition of substantial amounts of missing energy), yielding final 
states of the form xx + yy (+ , where in this work x and y are taken to be distinct 
SM fermions. Given a sufficient number of events, the observed topology (a 'box' or 
a 'wedge') clearly indicates whether or not Xi and Xj are identical particles. When 
simultaneous production of more than one pair of new particles is possible within 
a model, a more extensive set of topologies constructed from boxes and wedges 
(possibly with overlayed 'stripes') is obtained. A likelihood function indicating how 
well the a set of data points fits each possibility can be readily constructed if visual 
inspection does not suffice. The particular set of shapes the data sample should be 
thus compared to is of course model dependent. 

Though we wish to stress the general applicability of this technique to a fairly 
wide range of beyond-the-SM scenarios, application to i?-parity conserving SUSY 
models readily springs to mind. Thus the pair production of heavy MSSM neutralinos 



(excluding the lightest one, the LSP), with the subsequent decay of each -ino into 
a pair of leptons to aid identification, has been examined in detail. Here a fairly 
sizable number of distinct topological shapes is obtainable. This work then further 
specialises to -ino pairs produced in gluino/squark decays, most likely to be the 
dominant mode of -ino production at the LHC — if gluinos and/or squarks are 
relatively light. The number of possible topologies may be substantially reduced when 
this is the production mechanism compared to the E W production mechanisms which 
contain an -ino-ino-S" vertex if squark production does not re-introduce complexity. 
This was examined in some detail including possible tests of simulation results that 
may indicate the significance of squark production (and distinguish it from gluino 
pair production). Neutralino results thus obtained might be compared to those from 
charge asymmetries possible in samples of like-sign dilepton events from chargino 
pair production [T2"] . 

The 'hard edges' seen in a Dalitz-like plot yield information on the -ino mass 
differences as well as the identities of -inos participating in the decays (though it 
should be emphasised that the endpoints certainly need not equal the mass differ- 
ence of two -inos if on-mass-shell sleptons are involved in the decay chains), while 
comparing the relative densities of regions populated by different X t Xj production 
channels or combinations of channels has the potential to provide information on the 
relative production cross sections times leptonic BRs of these channels. We found 
simulation results from HERWIG for three distinctive points in MSSM parameter 
space (including cuts that nearly eliminate the backgrounds and a realistic detec- 
tor simulation) clearly closely tract the partial results we obtained at these points by 
'hand'-calculations based on the ISASUSY inputs. It is apparent from the Dalitz-like 
plots shown that this includes a substantial amount of information not available from 
a one-dimensional plot that just lumps together e + e~ and fi + fi~ invariant masses. 
Upcoming consideration of -ino pair production via heavy MSSM Higgs bosons 6], 
which also can have quite substantial rates in favourable regions of MSSM parameter 
space, will further expand the extra information obtainable from the two-dimensional 
Dalitz-like plot and thus should prove very exciting. Of course such analyses only 
incorporate mixed leptonic decays 17 - (e + e~[i + fi~ events, but not e + e~e + e~ or 
H + fi~ jjL + events). 

We also note that it is possible to make lego-style 3-dimensional plots with M(xx) 
and M(yy) along two axes and the binned number of events along a third axis. 
Figures obtained in this way were not found particularly illuminating for the specific 
processes and MSSM parameter points studied here (and with the modest amount 
of integrated luminosity assumed), but may be more useful in other studies. 

How far this method can go toward aiding reconstruction of the -ino mass spec- 
trum will depend on particulars of the point in the MSSM parameter space nature 
chooses, but clearly very significant information may be extracted. Given that an 
e + e~ linear collider with a centre-of-mass energy beyond that of LEP 2 is not ex- 

17 Unless angular correlations between leptons can be exploited to say how four same-flavour 
leptons should be arranged into two pairs without prejudicing the distribution. Or one could just 
plot all possible opposite-sign pair combinations, such plots may at least be distinguishable for 
different -ino pairs. 



pected for some time, it is crucial to seek the optimal methods for disentangling the 
-inos produced at the LHC Further, information on the heavier -ino states may 
prove crucial in deciding the reach of a future linear collider to perform the more 
precise measurements surely required. 
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Appendix 

The schematic Dalitz-like plot shown in the right square of Figure El is a collection 
of 6 observables (there labeled as regions a, (3, ...k) from which the production times 
leptonic BR values for the various -ino pairs, = 2,3,4) may be extracted. 

First, a triangular distribution of events is assumed for each -ino -ino mode: 

Tij = fC j dx J dy xy (23) 

where /C is a normalization constant that will drop out of the calculation. Now each 
region of the Dalitz plot contains events attributable to one or more of the modes 
The six different regions therefore correspond to different combinations of the 
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and the 
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where a, b, c, ...s are numbers between and 1 which represent the fraction of events 
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-£■0,1,2 being the three kinematical endpoints in the figure. Elements in each column of 
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the matrix M. must sum to unity. Now defining x 
yields 

b = y 
f = l-x 
J = %y(z - y) 
n = (1 — x) 2 c 
• = 2(1 -z)(z-y) s = (l-zf 

The linear system of equations is now easily solved for the individual rates r^. 



a = x 
e = y 2 
i = 2y(l — x) 
m = 2y(l — z) 
q = (1 — x)(l — z) 



c = x 
9 = z-y 
k = l-z 
(1 -x)(z-y) 



(t) -d z = (f 

d = xy 
h = 2x(l — x) 
I = x(l — z) 
P = (z-y) 2 
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